Using moss as a bio-indicator to evaluate soil quality in litchi orchard

Moss, was frequently found growing in litchi orchards. However, less known about whether it can be used as a visual bio-indicator for evaluating soil fertility and health. Therefore, soil chemical and biological properties, microbial community structures and the metabolic functions of microbes in soils between moss- and non-moss-growth areas were analyzed using traditional and high-throughput sequencing technologies. The results showed that pH and the contents of and available phosphorus (AP) in moss growth areas were significantly lower than those of non-moss growth areas, but the contents of alkali-hydrolyzable nitrogen (AN) and available potassium (AK) were significantly increased. In comparison with the soil of the non-moss-growth area, the abundances of hypotrophic microorganisms, such as Chloroflexi, Cyanobacteria and WPS-2 enriched in the soil of the moss-grown area. Moreover, the proportions of eutrophic microorganisms, such as Proteobacteria and Firmicutes also declined in the soils of the moss-growth area. Furthermore, the metabolic pathways of soil bacteria and fungi were also degraded in the moss-growth area. All above results indicated that not only lower soil fertility, but also soil microbial diversity also declined in moss growth area which compared to those of non-moss growth area. In one word, moss can be considered using as a visual bio-indicator for representing soil degradation in litchi orchards.


Introduction
Litchi (Litchi chinensis Sonn.) is a tropical and subtropical fruit, and it is mainly distributed in latitudes between 23˚and 27˚, such as in southern China, northern Vietnam, and Malaysia [1]. At present, it is widely grown in more than 20 countries around the world, and it has become has become a popular fruit [2]. However, its yield fluctuates year after year, and it is difficult to achieve yield stability [3]. In severe cases, even zero production can be found [4]. Insufficient nutrition is one of the important reasons for low litchi yields [5]. Previous studies have found that the contents of nitrogen (N) and potassium(K) in soil always induce low litchi yields [6]. In addition, litchi development requires abundant phosphorus (P), K, and calcium (Ca) during late fruit development [7]. As is well known, K plays an important role in fruit development, quality formation, and stress and disease resistances [8][9][10] available P in soil is positively correlated with the amount of P in litchi leaves, which affects the litchi yield. An inappropriate ratio of N to K results in fewer flower spikes and severe flower and fruit drop; soil fertility is an important factor in regulating litchi yields [11,12]. However, soil fertility cannot be immediately assessed by producers in practice. In particular, it is difficult to evaluate soil quality during the production of fruits, including litchi, citrus, and most other fruits. Recently, we found that moss could grow on the soil surface like a carpet in a litchi orchard (Fig 1A and 1B).
Previous studies have confirmed that moss is an important ground covering for degraded or damaged ecosystems, and it plays a vital role in their ecological properties [13]. Moss affects soil properties; in particular, it alters soil fertility through the accumulation of total N, K or alkali-soluble N, available K, P, and organic matter in the soil [14,15]. For example, moss enriches the P in soil during moss growth [16]; i.e., the more P is absorbed by mosses, the less P is left in the soil [17]. Meanwhile, moss is very sensitive to its environment and easily alters the water content and the nutrients in the environment with its growth [18][19][20][21]. Previous studies have considered moss as a bio-indicator for soil and water conservation, vegetation changes, environmental monitoring, and forest integrity [22,23].
In addition, soil microorganisms are important parts of ecosystems; not only they are closely linked to the energy cycle of terrestrial ecosystems, but they also participate in the cycling of material above and below the ecosystem floor, leading to the stability of the terrestrial ecosystem [24]. Bacteria are the most abundant taxa in soil, accounting for 70-90% of all soil microorganisms [25], and they play an important role in nutrient cycling and in maintaining soil health [26]. Fungi are also key members of soil microbiota; they are a major source of soil microbial biomass and are important drivers in maintaining cycles of matter in ecology [27]. The soil microbial community structure reflects the evolution of the ecosystem function, which relates to soil properties and environmental factors [9,28]. It has been confirmed that changes in soil properties often lead to alterations in the microbial community structure [29,30]. Moreover, soil microorganisms are very sensitive to the changes in their living environment, including soil nutrients [31], pH [32], soil physical properties [33], and plant diversity, among which soil nutrients and water content are the two most important factors affecting the diversity and function of soil microorganisms [34]. Additionally, soil microbial biomass C, N, and P are important indicators of soil health [35], and they can be used to evaluate soil quality, and crop productivity [36]. Soil enzyme activity is also an important index for soil quality assessment [37], and it can reflect the energy metabolism and material cycle in soil. As is wellknown, the soil microbial community structure is sensitive to environmental changes, and also can be used as a bio-indicator to reflect soil nutrient status and health. However, the bio-indicators of soil enzymes and soil microbial community structure, etc., are need to be analyzed firstly in lab, they are all not the visual indicators for soil quality.
Till now, studies on the relationship between moss and microbial community have focused on bacteria [38], Endophytic bacteria [39]. But minimal information is available on moss whether can be used as a visual bio-indicator for reflecting soil quality. Therefore, our aim here is to elucidate how differences of the soil biological properties, including soil enzymes, soil microbial biomass and the soil microbial compositions and their metabolic functions between moss and non-moss growth areas in a litchi orchard.

Field site description
The experimental site is located in a litchi orchard in Gao Li Village, Qinzhou City, Guangxi Zhuang Autonomous Region (109˚16 0 31@ E, 22˚12 0 30@ N), China. This area has a hilly topography and a subtropical monsoon climate; its average annual precipitation, temperature and frost-free period are 1658 mm, 22˚C and 348 d, respectively. Min and Max temperatures are 0˚C and 32˚C, respectively.

Soil sampling
The soil samples were collected in a Litchi orchard at Gaoli village (E: 109˚1621, N: 22˚1230), Qinzhou city of Guangxi Province, China on 26 May 2018. Feizixiao, a widely planted litchi cultivar, with 12-years old trees in the orchards under conventional managements were used in this study. And each litchi tree was applied 0.6 kg nitrogen, 0.18 kg phosphorus and 0.6 kg potassium using urea (N: 46.0%), superphosphate (P 2 O 5 :12.0%) and potassium chloride (K 2 O: 60.0%), respectively. Impurities on the soil surface were lightly removed, and soil samples from moss (hereafter abbreviated as M) and non-moss (hereafter abbreviated as L) areas were both carefully collected using a sterilized auger at a depth of 20 cm. Then, the soil samples were put into sterilized plastic bags and placed in an ice box with ice bags. After being transferred to the laboratory, soil samples were sieved through a 2 mm mesh stainless steel sieve, and then they were stored in a refrigerator at 4˚C for immediate analysis or were stored at −80˚C for later use. Meanwhile, parts of the soil samples were air-dried at room temperature for soil chemical analyses.

Soil chemical and biological properties
First, soil pH was measured using a pH meter (soil-water ratio 1:2.5). Total nitrogen (TN) was determined using the semimicro-Kjeldahl method [40]. Total phosphorus (TP) was determined by the alkali fusion-molybdenum anti-colorimetric method [41]. Total potassium (TK) was determined by alkali fusion-flame spectrophotometry [41]. Available P was determined using acid-fluoride solutions method [42], and available N, K were determined alkali diffusion method and flame photometry respectively [41].

Soil microbial biomass
The soil microbial biomass carbon (MBC) were determined using the 0.5 M K 2 S0 4 chloroform fumigation-extraction method as described by Vance et al. [43]. The soil microbial biomass nitrogen (MBN) was measured using the Ninhydrin-reactive nitrogen measurements of microbial biomass in 0 5 M K 2 SO 4 soil extracts method described by Joergensen et al. [44]. The soil microbial biomass phosphorus (MBP) was determined using fumigation-extraction with 0.5 M NaHCO3 at pH 8.5 as described by Powlson et al [36].

Soil enzyme activities soil microbial biomass
β-Glucosidase activity was determined using colorimetric estimation of ρ-nitrophenol released by β-glucosidase activity when soil is incubated in McIlvaine buffer (pH 4.8) with ρnitrophenyl β-glucoside and toluene at 30˚C for 1 hr method as described by Hayano et al. [45]. Aminopeptidase activity was measured using the using benzyloxycarbonyl phenylalanyl leucine as substrate, separated from the proteases by precipitation with 0.1 m CaCl 2 . And using Trisborate buffer extracts of a highly-organic method as described by Ladd et al. [46]. Acid phosphatase activity was determined using the universal buffer of pH 6.5, which containing 50 mM p-nitrophenyl phosphateat 37˚C for 1 h. The potential enzyme activities were defined as μg of p-nitrophenol (pNP) released by per g of soil within 1 h chloroform fumigation-extraction method as described by Tabatabai et al. [47].

Analysis of soil microbial community structure
Microbial community genomic DNA was extracted from samples using an E.Z.N.A.1 soil DNA Kit (Omega Bio-tek, Norcross, GA, USA) according to the manufacturer's instructions. The DNA extract was checked on 1% agarose gel, and the DNA concentrations and purity were determined using a NanoDrop 2000 UV-vis spectrophotometer (Thermo Scientific, Wilmington, NC, USA). PCR amplification and sequencing of the total DNA extracted from the rhizosphere soil samples were performed by Shanghai Majorbio Bio-pharm Technology Co., Ltd. (Shanghai, China), while PCR amplification was performed using an ABI GeneAmp1 type 9700 (ABI, CA, USA) and the products were recovered by 2% agarose gel electrophoresis, purified by an Axy PrepDNA gel recovery kit (AXYGEN) (Axygen Biosciences, Union City, CA, USA), eluted by Tris HCl, and quantified by the QuantiFluor™-ST (Promega, USA). Blue Fluorescence quantitative system (Promega). The purified amplicons were pooled in equimolar quantities and were paired-end sequenced (2 × 300) on the Illumina MiSeq platform (Illumina, San Diego, CA, USA) according to the standard protocols of Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China). Raw reads were deposited in the NCBI Sequence Read Archive (SRA) database (Table 1).

PICRUSt gene function prediction
PICRUSt was used to remove the effect of the number of copies of the 16S marker gene in the genome of the species and to standardize the OTU (Operational Taxonomic Unit) abundance table, using the green gene ID corresponding to each OTU. The KEGG Orthology (KO) information and COG family information corresponding to each OTU were obtained; then, the abundance of each COG and KO could be calculated. The COG database was parsed against Table 1. Sequencing type and primer sequence.

Source
Primer Name Primer Sequence Sequencing Platform Soil bacteria 338F 5 0 -ACTCCTACGGGAGGCAGCAG-3 0 PE300 the eggNOG database to obtain the descriptive information of each COG, as well as its functional information.

Statistical analyses
The trial data were statistically analyzed using Excel 2019 and Statistical Product and Service Solutions (SPSS) Statistics 21, and the results are shown as means with their standard deviations (mean ± SD). An online data analysis was performed using the free online Majorbio Cloud Platform (http://www.majorbio.com) of Majorbio Bio-Pharm Technology Co., Ltd.
(Shanghai, China). The variance inflation factors (VIFs) of different sample environmental variables were calculated, and redundant variables were removed. The significance was based on 999 Monte Carlo permutations. A linear discriminant analysis (LDA) and a linear discriminant analysis effect size (LEfSe) method were used to identify significantly different bacterial and fungal communities under different environmental samples [48].

Soil physicochemical properties
As shown in Table 2, the soil AP content and pH in the non-moss-growth area were significantly higher than those of the moss-growth area. On the contrary, the contents of soil AN and AK in moss-growth area were significantly higher than those of the non-moss-growth area. Meanwhile, the contents of soil TN, TP, TK and SWC were not significant difference between moss and non-moss growth areas (p > 0.05) ( Table 2).

Soil microbial biomass
As seen in Table 3, only MBN in the moss-growth area was significantly higher than that of non-moss growth area. In contrast, MBC and MBP were not significantly different between moss and non-moss growth areas (Table 3).

Soil enzyme activities
β-Glucosidase, aminopeptidase, and phosphatase are closely related to soil C, N, and P cycles, respectively. The results of the enzyme activities in soils of the moss and non-moss-growth Note: All data are presented as means ± standard deviation (SD). Different letters in the same column indicate significant differences among treatments at p < 0.05. SWC, soil water content; TN, total nitrogen content; TP, total P content; TK, total potassium content; AN, available nitrogen content; AP, available phosphorus content; AK, available potassium content. The same abbreviations are used in the tables below.
https://doi.org/10.1371/journal.pone.0278303.t002 Table 3. Soil microbial biomass C, N, and P between moss-(M) and non-moss growth(L) areas in litchi orchard (mg�kg −1 ). areas in litchi orchard were as follows. The activity of β-glucosidase was significantly higher in soils of the moss-growth area than those of the non-moss-growth area. By contrast, the activity of aminopeptidase was significantly higher in soils of the non-moss-growth area than those of the moss-growth areas However, the activities acid phosphatase between moss and non-moss growth area were not significant differences between each other (Table 4).

Correlation between soil bacterial community and the environmental factors
Based on the variance inflation factor (VIF) analysis, the environmental factors with p > 0.05 or VIF > 20 were screened and removed. First, as the VIF values of available potassium (AK) and total phosphorus (TP) were greater than 20. Therefore, they were removed. And the remained environmental factors were used for a correlation heatmap analysis. The results showed that Proteobacteria and Acidobacteria were positively correlated with available phosphorus (AP) in soil, which indicated that higher contents of AP are suitable for their growth.
Meanwhile, there was also a significant positive correlation between Cyanobacteria with the contents of total nitrogen (TN) and available nitrogen (AN). Moreover, the content of the total potassium (TK), it was positively and negatively related to Chloroflexi and Actinobacteria, respectively. Furthermore, Proteobacteria was positively correlated with pH And Cyanobacteria was positively correlated with SWC. A linear regression analysis also showed that TN (R 2 = 0.8231, p = 0.0125), TK (R 2 = 0.3, p = 0.2606), AP (R 2 = 0.0162, p = 0.8099), AN (R 2 = 0.1331, p = 0.4771), pH (R 2 = 0.1498, p = 0.4485), and SWC (R 2 = 0.2445, p = 0.3074) which suggested that soil bacterial community structure could be altered by the main environmental factors, particularly, the content of the total nitrogen was the most significant factor among the above environmental factors (Fig 2).

Alpha and Beta diversities of soil bacterial community
The Shannon and Simpson indices (Fig 3A and 3B), which describe soil bacterial diversity, did not show significant differences between moss and non-moss growth areas. Moreover, Ace and Chao1, the bacteria richness indices, also showed the same trends between moss and nonmoss growth areas. These results suggested that soil bacterial diversity and richness were not significantly changed in moss growth areas which compared to the non-moss growth areas in litchi orchard (Fig 3C and 3D). Based on the Bray-Curtis distance from PCoA, the bacterial community was significantly aggregated into two groups corresponding to the compositions of the soil bacterial communities between non-moss and moss-growth areas in litchi orchard (ANOSIM: R = 0.3333, p = 0.1990; Adonis: R 2 = 0.3336, p = 0.2000). At the OTU level, the contribution rates of the first principal component (PC1) and the second principal component (PC2) to the samples were 44.12% and 37.61%, respectively. Meanwhile, soil bacteria in the non-moss growth area grouped together which precisely separated with those of the moss growth area in litchi orchard (Fig 4). Table 4. Soil enzyme activities related to C, N, and P cycles between moss (M) and non-moss-growth (L) areas in litchi orchard (nmol�g −1 �min −1 (30˚C)).
All above results suggested that soil bacterial compositions at the phylum level in moss growth areas were significantly different with those of the non-moss growth area in litchi orchard (Fig 5A). At the genus level, the numbers of soil-dominant bacterial genera (i.e., relative abundances greater than 1%) in the moss and non-moss growth areas in the litchi orchard were 22 and 17, respectively. The dominant soil bacterial genera in the moss-growth area were norank_c_AD3 (12.52%), norank_o_Chloroplast (7.01%), norank_f_Xanthobacteraceae (6.24%), Acidothermus  (Fig 5B).
At the genus level, 412 and 441 soil bacterial genera could be detected between the moss and non-moss-growth areas in the litchi orchard, respectively. Among them, there were 353 common soil bacterial genera between the moss and non-moss-growth areas. In addition, there were 59 and 88 unique soil bacterial genera in the moss-growth and non-moss-growth areas, respectively ( Fig 6A). Moreover, at the OTU level, there were 2214 and 2376 soil bacterial genera in the moss and non-moss growth areas in the litchi orchard, respectively. Among them, 1727 common soil bacterial genera were found in the moss and non-moss growth areas. Moreover, there were 649 unique bacterial OTUs in the moss growth area and 487 unique bacterial OTUs in the non-moss-growth area (Fig 6B).
Significant differences of soil bacteria between the moss and non-moss growth areas in litchi orchard and the main contributing biomarker classes were examined using an LEfSe analysis (LDA threshold of 3.5) (Fig 7). A total of 51 soil bacteria clades showed statistically significant differences between moss and non-moss growth areas. Solirubrobacterales (from order to family), Leptolyngbyales (from order to family), Ktedonobacteria (from class to genus), norank_f__Ktedonobacteraceae (genus), 1921-2 (genus) and WPS-2 (from phylum to genus) enriched in the moss growth area.

PLOS ONE
Using moss as a bio-indicator to evaluate soil quality in litchi orchard

Functional predictions of soil bacterial community structures
Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, the primary functional layer of soil bacteria in the moss and non-moss-growth areas contained six types of biometabolic pathways (Fig 8), i.e., Metabolism, Environmental Information Processing, Cellular Processes, Genetic Information Processing, Human Diseases, and Organismal Systems.
As seen in Fig 8, the percentages of the bio-functional pathways in the primary layer of the soil bacteria in the moss-growth area were inferior to those of the non-moss-growth area.
In addition, the secondary functional layer of soil bacterial genes in the moss and non-mossgrowth areas mainly consisted of Membrane transport; Cellular community-eukaryotes; Signaling molecules and interaction; Circulatory system; Drug resistance: Antimicrobial; Cell motility; Antimicrobial eukaryotes; Signaling molecules and interaction; Circulatory system, etc. (Fig 9).
In Cell motility, Neurodegenerative disease, Infectious disease, viral, Cancer, specific types, Infectious disease, parasitic, Circulatory system and Cellular community-eukaryotes and others 40 secondary functional layers of soil bacteria in the moss growth area were degraded. However, the functional predictive gene types of Excretory system, Digestive system, Immune disease, Development and regeneration, Sensory system, and Signaling molecules and interaction showed an increasing trend in the moss growth area. Particularly, in comparison with non-moss growth area, the Signaling molecules and interaction increased from 38.59% to 61.41% in the moss growth area. The above results indicated that although the predicted gene types in the secondary functional layer were similar between the moss and non-moss growth area, but the gene copy numbers were significant differences between each other. The functional degradation of soil bacteria in the moss growth area indicated that the functional abundances of soil bacterial genes at the secondary functional layer were lower than those of the non-moss growth area in the litchi orchard.

Correlation between soil fungal community and the environmental factors
The environmental factors with p > 0.05 or VIF > 20 were screened using a variance inflation factor (VIF) analysis. As the VIF values of the environmental factors, available potassium (AK) and total phosphorus (TP) were greater than 20, therefore, they were removed. Meanwhile,

Alpha and Beta diversities of soil fungal community
The Shannon and Simpson indices (Fig 11A and 11B) did not show significant differences in fungal diversities between moss and non-moss growth areas. Moreover, the Ace and Chao1 indices, also showed the same trends with diversity between the moss and non-moss growth areas (Fig 11C, 11D). These results indicated that soil fungal diversity and richness in moss growth area did not significantly changed in litchi orchard.
Based on the Bray-Curtis distance from PCoA, the soil fungal communities were significantly clustered into two groups between moss and non-moss growth areas in litchi orchard (ANOSIM: R = 0.4815, p = 0.195; Adonis: R 2 = 0.3304, p = 0.100). As seen at Fig 12, the contribution rates of PC1 and PC2, as the first and the second principal components were 44.48% and 20.05%, respectively (Fig 12).

PLOS ONE
and Mortierellomycota decreased, but the proportion of Ascomycota increased in moss growth area. Moreover, unclassified_k_Fungi, as the soil dominant fungal phylum in the nonmoss-growth area, but it lost in the moss growth area. Furthermore, Rozellomycota was the unique soil dominant fungal phylum of the moss growth area in litchi orchard (Fig 13A).
Meanwhile, unclassified_c_Leotiomycetes, Phoma, Geminibasidium, Tremella, Omphalina, Rickenella, and unclassified_o_GS11t were the unique soil dominant fungal genera in the moss growth area. In contrast, unclassified_o_Hypocreales, Fusarium, Trichoderma, unclassified_c_-Sordariomycetes, unclassified_k__Fungi, Apiotrichum and Acaulium were the special soil dominant fungal genera in non-moss growth area in litchi orchard. Moreover, 302 and 337 soil fungal genera were obtained in the moss and non-moss growth areas in the litchi orchard, respectively. Among them, 238 common soil fungal genera were found between moss and non-moss-growth areas. Furthermore, 64 and 99 unique fungal genera were also detected in the moss and non-moss growth areas, respectively (Fig 14A). At the OTU level, 1238 and 1484 soil fungal OTUs were found in moss and non-moss growth areas, respectively. Among them, 790 OTUs were common fungal OTUs, and 448 and 694 were unique fungal OTUs in the moss and non-moss growth areas, respectively ( Fig 14B). All above results suggested that soil fungal compositions had been changed in the moss-growth area. In comparison with the non-moss-growth area, not only at the genus level, but also at OTU levels, a degraded trend of the fungal compositions could be found in the moss-growth area in the litchi orchard.
Significant differences between moss and non-moss-growth areas in litchi orchard and the main contributing biomarker classes were examined using the LEfSe analysis (LDA threshold of 3.5) (Fig 15).
The compositions of the fungal communities differed significantly between the non-moss and moss-growth areas in litchi orchard according to the LEfSe analysis. Glomerellales (from order to genus), Leotiomycetes (from class to genus), Thyridariaceae (from family to genus), Didymellaceae (from family to genus), Hymenochaetales (from order to genus), and Pucciniomycetes (from class to genus) enriched in soil of the moss-growth area. In contrast, Sordariomycetes (class), Nectriaceae (from family to genus), Microascales (from order to genus), and Onygenales (from order to genus) enriched in soil of the non-moss-growth area in litchi orchard.

Functional predictions of soil fungal community structures
By comparing with the EggNOG database, 16 COG functions were obtained. The analysis of interorganizational COG function composition showed that the top 10 COG functions in the

Soil fertility indexes between moss and non-moss growth areas in litchi orchard
Soil enzymes involved in the cycling of all nutrient elements, their activities related to the C, N, and P cycles in soil. And they had been regarded as the important biological indicators for soil fertility and ecological diversity, stability assessment [49][50][51]. Meanwhile, soil microbial biomass also can be used as the biological indicator for reflecting energy cycling, nutrient transport in soil [52,53].
First, even though the contents of total nitrogen, phosphorus and potassium were not significantly different in soils between moss and non-moss growth areas. However, the contents of available nitrogen, phosphorus and potassium in soil were all significantly different between each other. Although the contents of available nitrogen and potassium in soil of moss growth area were significantly higher than those of non-growth area, but the content of available phosphorus significantly decreased in moss growth area. Meanwhile, in comparison with nonmoss growth area, soil pH in moss growth area was also significantly declined.
In addition, MBC and MBP were not significantly different between moss and non-moss growth areas, but the MBN in moss growth area was significantly higher than that of nonmoss-growth area. Moreover, except of the activity of acid phosphatase, the activities of β-Glucosidase and aminopeptidase in soils were all significantly different between moss and nonmoss growth area. i.e., the activities of β-Glucosidase and aminopeptidase in soil of moss growth area were significantly higher and lower than those of non-moss growth area, respectively. All above results suggested that soil fertility in moss growth area exactly had already changed. Particularly, soil acidification and phosphorus deficiency were the soil primary features of the moss growth area in litchi orchard.

Soil microbial compositions and functions between moss and non-moss growth areas in litchi orchard
Soil microorganisms can be used as important indicators to evaluate soil fertility and represent its changes [54]. Numerous studies have confirmed that soil environmental factors have a significant influence on soil microbial community structure [55]. Soil nutrient status also can affect the functional diversity and abundance of microbial populations in soil ecosystems [56]. As proteobacteria and firmicutes are eutrophic bacteria, which grow rapidly under eutrophic conditions [57,58]. However, the proportion of proteobacteria in soils of the moss growth area was 18.18% only, but it was 37.58% in the non-moss growth area. Meanwhile, Firmicutes was not the soil dominant bacterial phylum in the moss growth area, but it was the soil dominant bacterial phylum in the non-moss-growth area. Moreover, as soil fertility can be significantly improved by Mortierellomycota, which it can promote the contents of soluble carbon, available phosphorous and phosphorus-related enzyme activities [45,59]. And we also found that the proportion of Mortierellomycota in the non-moss growth area was also higher than that of the moss growth area. In addition, Chloroflexi and Cyanobacteria also are considered as the malnourished bacteria [60,61]; WPS-2 generally lives in low soil fertility [62]; and Rozellomycota also shows strong adaptations to extreme soil environments, such as extreme pH levels and soil temperatures [63]. However, in comparison with non-moss growth area, Chloroflexi, Cyanobacteria, WPS-2 and Rozellomycota were all the soil dominant microbes in moss growth area. All above results indicated that higher nutrients status in non-moss growth area could be inferred than those of -moss growth area.
In addition, although soil bacteria and their metabolic pathways of the primary functional level between moss and non-moss growth areas were similar, but some bacterial gene functions, such as the bio-functional pathway, its proportion was lower in the moss growth area than that of the non-moss growth area. Furthermore, most of the fungal functional group's proportions in moss growth area were also found lower than those of the non-moss growth area. The results suggested that non only soil microbial compositions degraded, but also soil microbial functions also declined in moss growth areas which compared to non-moss growth areas in litchi orchard.

Conclusions
In comparison with the soils of the non-moss growth area, soil pH and the contents of AP were significantly decreased in the moss growth area. Meanwhile, the abundances of hypotrophic microorganisms, such as Chloroflexi, Cyanobacteria and WPS-2 enriched in the soils of the moss-grown area. Moreover, the proportions of eutrophic microorganisms, such as Proteobacteria and Firmicutes also declined and the metabolic pathways of soil microbial functional groups were also degraded in the moss growth area in litchi orchard. In one word, moss can be considered using as a visual bio-indicator for representing soil degradation in litchi orchards.